Thermal fluctuations and flux-tunable barrier in proximity Josephson junctions 
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The effect of thermal fluctuations in Josephson junctions is usually analysed using the 
Ambegaokar-Halperin (AH) theory in the context of thermal activation. "Enhanced" fluctuations, 
demonstrated by broadening of current- voltage characteristics, have previously been found for prox- 
imity Josephson junctions. Here we report measurements of micron-scale normal metal loops con- 
tacted with thin superconducting electrodes, where the unconventional loop geometry enables tuning 
of the junction barrier with applied flux; for some geometries, the barrier can be effectively elim- 
inated. Stronger fluctuations are observed when the flux threading the normal metal loop is near 
an odd half-integer flux quantum, and for devices with thinner superconducting electrodes. These 
findings suggest that the activation barrier, which is the Josephson coupling energy of the proximity 
junction, is different from that of conventional Josephson junctions. Simple one dimensional quasi- 
classical theory can predict the interference effect due to the loop structure, but the exact magnitude 
of the coupling energy cannot be computed without taking into account the details of the sample 
dimensions. In this way, the physics of this system is similar to the phase slipping process in thin 
superconducting wires. Besides shedding light on thermal fluctuations in proximity junctions, the 
findings here also demonstrate a new type of superconducting interference device with two normal 
branches sharing the same SN interface on both sides of the device, which has technical advantages 
for making symmetrical interference devices. 

PACS numbers: 74.45. +c, 74.50.+r, 74.40. -fk, 74.78.Na 



Understanding the effect of thermal fluctuations in 
nanoscale superconducting devices is important for 
applications.1' One particular type of nanoscale super- 
conducting device is a proximity junction,'212J consisting of 
a normal metal wire contacted by two superconducting 
wires as electrodes. The difference between such a de- 
vice and a conventional Josephson junction made of bulk 
superconducting electrodes separated by a thin insulator 
is twofold: 1) Instead of tunnelling through the insu- 
lating barrier, charge carriers diffuse through the normal 
metal barrier in proximity junctions, leading to a stronger 
coupling between the two superconducting electrodes for 
proximity junctions; 2) the thin superconducting elec- 
trodes of these nanoscale junctions are different from bulk 
superconductors since their dimensions are now compa- 
rable to the superconducting coherence length S,s- 

These differences may explain the previously observed 
strong fluctuations in nanoscale proximity junctions,'^"^ 
which manifest themselves as "enhanced" broadening of 
the current-voltage characteristics (CVC) compared to 
the "intrinsic" broadening of CVC due to thermal activa- 
tion of a resistively shunted junctio n ( RSJ), as described 
by Ambegaokar and Halperin (AH). EH! This enhancement 
can be characterized by an effective noise temperature 
Tjv, which is higher than the bath temperature T^^ 

In this paper, we report measurements of a special type 
of proximity junction with a loop structure embedded 
into the junction itself. By threading a magnetic flux 
through the loop, we find T/v is maximum when the flux 
through the loop is an odd half-integer of the supercon- 
ducting flux quantum $o = /i/2e, for fixed T;,. The ratio 



between T/v and T;, depends on the geometries of the 
particular sample, but does not depend on the bath tem- 
perature Tf,. The fiux dependence of T/v is better under- 
stood if we consider a phase slipping process similar to 
that in thin superconducting wires, which again suggests 
that for a nanoscale proximity junction, the activation 
energy is different from the standard Josephson coupling 
energy. Such devices exhibit many of the properties of 
dc SQUIDs, but with the advantage that the devices can 
be designed to be almost perfectly symmetric, as the two 
normal branches share the same SN interface on both 
sides of the device, hence allow unprecedented tunability 
of the system by means of an external magnetic flux. 

The remainder of this paper is organized as follows: In 
Section I, a general introduction for thermal activation 
in superconducting wires and Josephson junctions is pre- 
sented, followed by a discussion of how quasiclassical the- 
ory can be used to extend the thermal activation model 
to proximity junctions. Then the experimental results of 
magnetoresistance (Section II), the current voltage char- 
acteristics (CVC) (Section III), and temperature depen- 
dence of resistance (Section IV) are discussed. Finally 
we summarize the flndings and propose further theoreti- 
cal study beyond one dimensional quasiclassical theory. 



I. THEORETICAL BACKGROUND 

The dynamics of a particle trapped in a shallow po- 
tential well is a fundamental problem that has wide ap- 
plicability to a number of areas in statistical physics. At 



finite temperatures, due to coupUng to a thermal bath, 
the particle may escape from the potential well through 
the process of thermal activation over the barrier repre- 
sented by the edge of the potential. 

For superconductors, two specific phenomena have 
been explored extensively in terms of the physics of ther- 
mal activation. The first phenomenon is the generation 
of phase-slips in a thin superconducting wire with cross- 
section siz e com parable to the superconducting coherence 
length ^g' ^oi^^l In this case, thermally activated phase 
slips (TAPS) lead to the appearance of a resistance tail of 
the thin wire at temperatures below the nominal transi- 
tion temperature Tc. The voltage V generated by TAPS 
is related to the time evolution of the macroscopic su- 
perconducting phase (p through the Josephson relation 
2eV/h — dip/dt, which on average is determined by the 
number of TAPS per unit time, each TAPS correspond- 
ing to a change of 27r. At the instant in space and time 
where such a phase slip event occurs, the superconduct- 
ing order parameter vanishes, which costs an energy A.F, 
the barrier over which the system must be thermally ac- 
tivated. In the earlier pioneering theory by Little,^° the 
relevant energy barrier is the condensation energy in a 
small coherent volume determined by the dimensions of 
the superconducting wire, and the temperature depen- 
dence of resistance R(T) is described by an Arrhenius 
type equation. Later developments gave a more accurate 
description of AF and the attempt frequency for acti- 
vation oyer the barrier, often referred to as the LAMH 
theory.lSIIll 

The second phenomenon is the onset of finite voltage 
in a Josephson junction. A resistively shunted junction 
(RSJ) can be modelled as a particle in a one-dimensional 
washboard potential in a viscous medium, where the 
distance coordinate corresponds to the phase difference 
if across the junc tion, as shown by Ambegaokar and 
Halperin (AH) ^^ With no current through the junction, 
the system sits in a local minimum of the potential. Ap- 
plication of a current through the junction corresponds 
to tilting the washboard, shifting the position of the local 
minima. At some value of current less than the nominal 
critical current Ic, the washboard potential is tilted suffi- 
ciently for the system to be thermally activated over the 
barrier between two adjacent potential minima. Once 
this occurs, the system will continue to roll down the 
washboard potential, corresponding to a continuous time 
evolution of (p, and hence a finite voltage will appear 
across the junction according to the Josephson relation. 
Here the energy barrier (the height of the washboard po- 
tential) being thermally activated over is the Josephson 
coupling energy Ej = {h/2e)Ic, where Ic is the critical 
current. 

A proximity junction made by a normal metal wire be- 
tween two superconductors exhibits properties different 
from a conventional tunnelling junction. Proximity SNS 
junctions have also been investigated for a long time.-'^ 
More recent theoretical investigations have used the 
framework of the quasiclassical theory of superconduc- 



tivity to discuss the characteristics of SNS junctions.'Sl 
Qualitatively, the physics of the SNS junction is well 
understood. The proximity to the superconductor has 
two major effects on the quasiparticles in the normal 
metaLliS First, it induces superconducting-like correla- 
tions between quasiparticles that increase the conduc- 
tance of the normal metal, and second, it induces a gap 
in the quasiparticle density of states A^(£').'i^ While in 
the superconductor the energy scale for the density of 
states is given by A, the energy scale for N{E) in the 
normal metal for the diffusive case is set by the Thou- 
less energy Erh = HD/L'^, where D = vpi/S is the 
quasiparticle diffusion constant in the normal metal, £ 
being the elastic scattering length, and L is the length of 
the normal metal. Similarly, the maximal supercurrent 
that can flow through such a proximity junction is set 
by Eth, not by A as for a conventional Josephson tun- 
nel junction, when Exh ^ A (the long junction limit). 1^ 
However, it is difficult to obtain quantitative predictions 
for the current-voltage characteristics except in limiting 
cases. The problem arises from the fact that a finite volt- 
age results in a time dependent phase difference between 
the superconducting electrodes, which in turn leads to 
time-dependent boundary conditions for the quasiclas- 
sical equations. These equations have been solved un- 
der certain simplifying assumptions, such as low inter- 
face NS interface transparency,!^ or in the limit of volt- 
ages small compared to Eth/s!^ Unfortunately, these 
assumptions do not apply to our samples. Consequently, 
in the discussion below, we have chosen to discuss our re- 
sults phenomenologically in the framework of RSJ model 
discussed above, but with the energy potential landscape 
being calculated using the quasiclassical theory of super- 
conductivity. 

The Josephson junction model leads to an interesting 
extension, where two junctions can be connected in par- 
allel to form a dc superconducting quantum interference 
device (SQUID). In this case, there are two independent 
parameters (the phase differences across the two junc- 
tions), giving rise to a two-dimensional potential for the 
system. The system can transition from one local mini- 
mum to another through saddle-points in the potential.'^Sl 
Due to the fact that the position of the system on the two 
dimensional potential is sensitive to the external mag- 
netic flux <f>, dc SQUIDs have bee n inve stigated exten- 
sively due to their device potential. ^-'^ ^^ With any real 
dc SQUID, the two Josephson junctions cannot be fabri- 
cated to be exactly the same, which restricts the ability 
to tune the system. 

As with conventional junctions, two proximity junc- 
tions may be combined in parallel to form a dc SQUID. 
Schematics of two types of SNS junctions are shown in 
Fig. [l] In the asymmetric device shown in Fig. 1(a), the 
normal metal arms (shown in gold) connect to the su- 
perconducting wires (shown in gray) at different points, 
likely resulting in different NS interface transparencies for 
the two arms of t he d evice, and consequently an asym- 
metric dc SQUID .EESl However, the long-range nature of 




°'5 - 1 

*/*„ -0.5 -2 '^ 



*/*„ 0,5 -2 "^ 





0/it 



FIG. 1: (Color online) Schematic diagrams of two SNS quan- 
tum interference devices: asymmetric (a) and symmetric (b). 
The normal metal arms are shown in gold and the supercon- 
ducting wires are shown in gray. The arrow in both figures 
corresponds to the the direction of applied magnetic flux. In- 
sets: Scanning electron micrographs of the devices measured, 
the scale bars are 1 fim. The calculated energy profile of the 
asymmetric/symmetric device are shown in (c)/(d) as a func- 
tion of the external magnetic flux $ and the phase difference 
across the two superconductors (j>, which is determined by the 
external current through the device. The energy profiles are 
calculated based on the quasiclassical theory, as described in 
the text. 



the Josephson coupling in SNS devices enables a new 
type of device that is not possible with conventional tun- 
nel junctions. As shown schematically in Fig. fl|b), the 
device consists of a single normal metal loop between 
the two superconducting contacts, so that the NS inter- 
face transparencies for the two arms of the device are 
the same. In this device, the modulation of quantum 
interference by an external magnetic flux occurs within 
the junction itself, i.e., the superconducting phase wind- 
ing happens along the loop inside the junction, similar 
to the case of a superconducting loop that shows clas- 
sical Little-Parks oscillations .I^St^ll As the NS interfaces 
are the same for both arms of the loop, this device po- 
tentially can behave as a perfectly symmetric dc SQUID, 
provided that the length of the two normal metal arms 
are the same. 

To understand the flux-tunable thermal activation bar- 
rier in these SNS junctions, we use quasiclassical theory 
and a simple one dimensional model to calculate the en- 
ergy of the system as a function of the phase difference ip 
between the two superconductors (see appendix for de- 
tails). The energy of the system is given by (see, e.g., 
Ref. [m, page 198) 



Eji^,^) 



h 

2e 



V/,(<^',$), 



(1) 



where /^((/j',^) is the supercurrent through the system, 
which is a periodic function of the phase difference ip' 



and the externally applied flux $. To calculate Is^ we 
use the extended circuit theory^SI ^j^,^ numerically solve 
the Usadel equations for the sample geometries shown in 
Figs. llTa) andllTb). In the simple model we assume that 
the interfaces between the normal metal (N) and the su- 
perconductors (S) are perfectly transparent, the gap A 
regains its bulk value in the superconductor within a very 
short distance of the NS interface, and that the distribu- 
tion of quasiparticles in the superconducting reservoirs 
is given by the equilibrium Fermi function f{E). Since 
the characteristic unit for the supercurrent for an SNS 
junction is Erh/eR, we use {hl2e'^R)ETh as the charac- 
teristic unit of energy for EJ^^^M\t}l L being the length 
of one side of the loop. For the parameters used in this 
simulation, the amplitude of the supercurrent is about 
0.2 Exh/eR, so the modulation of Ej at fixed $ is about 
0.4 {h/2e^R)ETh (note that since Ej has an arbitrary 
constant from the integration in Eq. (1), in Fig. n] we 
assume Ej = a.t ip = — 27r). 

The resulting energy profiles are shown in Figs, [lie) 
andllrd) as a function of ip and $, the two parameters un- 
der external experimental control. If $ is fixed at integral 
values of the superconducting flux quantum $o = h/2e, 
there is an energy barrier for evolution of the phase, as 
shown by the trajectories of the blue particles at zero 
fiux, so that at low temperatures, the phase ip is station- 
ary, and no voltage is developed across the device. For 
half-integral values of the applied fiux (^ = {n + l/2)$oi 
where n is an integer) there is a difference between the 
asymmetric and symmetric cases. For the asymmetric 
case, there is still a small energy barrier, as shown by 
the trajectory of the red particle in Fig. flTc). Conse- 
quently, the resistance at odd half-integral flux quanta 
will eventually vanish if the temperature is low enough, 
in the absence of quantum tunnelling. In contrast, for 
the symmetric case, there is no energy barrier at odd 
half-integral flux quanta (as shown by the trajectory of 
the red particle in FigflTd)), so that the device will have 
a finite resistance even at the lowest temperatures. 



II. ZERO BIAS MAGNETORESISTANCE 

Figures l2ja) andl2jb) show the resistance of the asym- 
metric and symmetric devices respectively as a function 
of applied magnetic fiux at a number of temperatures. 
The details of fabrication and measurement are similar 
to those reported elsewhere.'^l The data are taken in the 
limit of zero dc current, with only a very small ac current 
for the resistance measurement (about 10-20 nA). 

At higher temperatures, the resistance is finite at all 
values of $ for both geometries, and is periodic in the 
applied flux, with a fundamental period of $0- As the 
temperature is lowered, the resistance for both devices 
vanishes near $ = n$o- As the temperature is lowered 
still further, the resistance of the asymmetric device also 
vanishes at odd half-integral values ^ — {n + l/2)$o of 
the applied flux. In contrast, the resistance of the sym- 
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FIG. 2: (Color online) Resistance as a function of applied 
magnetic flux for the asymmetric sample (a) and the symmet- 
ric sample (b) at 0.6 K (blue curves), 0.4 K (green curves), 
and 0.03 K (red curves). While the oscillations for the asym- 
metric device die out rapidly with decreasing temperature, 
the oscillations for the symmetric device survive to the low- 
est measurement temperature, (c) Magnetoresistance of the 
symmetric device around $ = $o/2 at four different temper- 
atures. The solid lines are fits to the AH theory as described 
in the text. Below 0.15 K, the magnetoresistance does not 
change with temperature, (d) Critical current at zero applied 
flux for the symmetric and asymmetric devices. The plus 
symbols show the critical current expected from the fits of 
(c), multiplied by a factor of 5, as described in the text. 



metric device at $ = (n + 1/2) $o remains finite down to 
the lowest temperatures, while the peak width narrows as 
the temperature is lowered, saturating below about 0.15 
K, as shown for the peak around $ = $o/2 in Fig. [2jc). 
According to the physical picture for the symmetric 
device in Fig. W[d) , there is a small but finite barrier for 
the system to overcome, except at exactly half-integral 
values of flux. This barrier decreases monotonically as 
the system approaches $ = (n -f l/2)<I>o- In fact, for 
conventional junctions the Josephson coupling energy is 
Ej — {h/2e)Ic, Ic being the critical current. In our SNS 
junction case, Ej is given by Eq. (flj), and is a function 
of the external flux $. Assuming symmetrical long junc- 
tions, we obtain the usual sinusoidal dependence of the 
current J, on ip, with the Josephson coupling energy 



E, 



2e ' 



(2) 



where Ic{^) = ^c(0)| cos(7r$/$o)|j the same as that for 
a symmetric dc SQUID (see, e.g., Ref. [23, page 215). 
Thus, as the temperature is lowered, the system needs 
to be closer to <I> = (n -|- 1/2) $o until the energy barrier 
is low enough for the particle to jump over and for the 
junction to exhibit a finite resistance. At even lower tem- 
peratures, the system may tunnel through the barrier. 



resulting in a temperature independent resistance. How- 
ever, it is necessary to verify whether the electron tem- 
perature follows the bath temperature or not ,1221 which is 
non-trivial and will not be discussed here. We note that 
the measured critical current saturates at lower temper- 
atures and is about 10 times smalle r tha n the value pre- 
dicted in the zero temperature limit ii2iS3i probably due to 
an imperfect SN interface, fluctuations, and heating. ^^ 

To quantitatively model the magnetoresistance around 
half-integral flux values for the symmetric device, we can 
use AH theory for RSJ model, as extended here to a dc- 
SQUID with a flux-tunable barrier. The AH theory pre- 
dicts that the normalized resistance in the zero current 
limit is given by 



RAH=Io'il/'^) 



(3) 



where Iq is the modifled Bessel function, and 7 = 
2Ej/kBT is the ratio between the barrier and the ther- 
mal fluctuation energy. From Eq. (pi), we have 7 — 
hlc{^) / eksT , and we can then flt the magnetoresistance 
of the symmetric device near <I> = $o/2 using the AH 
theory with Eq. ([3|. 

The solid lines in Fig. I2jc) show the resulting flts of 
R{^) at four different bath temperatures, using only the 
measured peak resistance -Rp at $ = $0/2 and zero fleld 
critical current /c(0) as fltting parameters. Figure [2jd) 
shows a comparison of /c(0) obtained from the flts (mul- 
tiplied by a factor of 5) compared to the experimentally 
measured values of /c(0). The fltted values of /c(0) are 
a factor of 5 smaller than the experimentally measured 
values of lc(0) over the entire temperature range. Con- 
versely, we can claim that the effective noise temperature 
T/v is 5 times larger than Tf, since only the ratio between 
Ej and ksT matters. 

The strong correlation between the measured and the 
fltted values indicates that the AH theory is somewhat 
applicable for the flux-tunable proximity junction near 
$ — <I'o/2, but the effective potential barrier for thermal 
activation appears to be smaller by about a factor of 5 
compared to Eq. (pi). As this factor is temperature in- 
dependent, this "enhanced" fluctuation could not be due 
to quantum fluctuations,^^ nor due to heating. In fact, if 
there is heating during measurements of Ic then the mea- 
sured Ic should be smaller than that inferred from the 
thermal activation model at zero current limit, but the 
opposite is observed. This discrepancy could be related 
to the assumption of a standard external shunt resistor as 
the thermal noise source in the RSJ model, while the SNS 
junction is self-shunted, or it could be due to the neglect 
of inverse proximity effect in superconducting electrodes 
and formation of minigap in the normal metal,'^ which 
may change the relation between Ic and Ej in Eq. (|2|. 
To better understand this "enhanced" fluctuation, below 
we characterize devices beyond the zero bias limit. 




1), the normalized voltage v — V/Ic'R is 
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FIG. 3: (Color online) Differential resistance at 0.4 K for si, 
(a), and for s2, (c). Different colors correspond to different 
applied magnetic flux: (magenta), 0.25 (blue), 0.375 (red), 
and 0.5 (green) in units of $o. Solid lines are fits using AH 
model, (b) and (d) show the flux dependence of the noise 
temperature Tn (red) and critical current Ic (blue) extracted 
from the AH flts of the differential resistance in (a) and (c). 
The insets are SEM images of the devices with false color 
enhancement. The area colored in gold is Au and the area 
colored in purple is Al. The white scale bars are 1 /im. 



TABLE I: Sample parameters of two representative proximity 
junction devices: wai and Iai are the width and thickness of 
the AI electrodes, wau and tAu are the width and thickness 
of the Au loops, Tc and He are the critical temperature and 
critical magnetic field of the Al electrodes, which is sensitive 
to Wai for quasi-lD wires. 
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III. FINITE BIAS 



Finite bias measurements were conducted on two rep- 
resentative proximity junction devices with supercon- 
ducting Al electrodes of different widths: about 110 nm 
and 220 nm respectively for devices si and s2, see Ta- 
ble [l] for other parameters, and insets of Fig. |3] for SEM 
images. The Al wires extend for several micrometers on 
each side of the loop before overlapping with Au leads 
for four-probe measurements. 

As in early investigations for conventional Josephson 
junctionSjI^^lEII-^ye measured the differential resistance of 
the two proximity junction devices at 0.4 K, as shown in 
Fig. |3] and fit the data by numerical differentiation of the 
voltage in the AH theory.^ In the limit of small currents 
{x = I/Ic < 1) and low temperatures (7 = Icfi/eksT 3> 



2v/r 



x-^ exp[— 7 



[/T 



' + X sin x]] sinh[7r7a;/2] . 

.(4) 
As for Rah in Eq. ([3|, the essential fitting parameter is 7, 
the ratio between the energy barrier Ej and the thermal 
energy kBTi,. Assuming Ej is just the Josephson coupling 
energy in Eq. ([2|, then we need to replace Ti, (0.4 K) 
with an effective noise temperature T/y to characterize 
the enhancement of fluctuations. 

In Fig. [3] we show the result of fitting at several dif- 
ferent flux values (in unit of $0), with both Tjy and Ic 
used as fitting parameters. When Ej is suppressed at 
around half flux quanta, Tpf reaches its maximum, and 
the quality of the fit is good to large values of /. When Ej 
is maximum at around integer flux quanta, there is clear 
deviation of the fit from the experimental data as the bias 
current approaches I^ and we can only fit the low current 
bias regime. Compared to previous investigationsj-^^^ 
here we can tune the Josephson coupling (Ej) without 
varying the physical temperature, which gives us a knob 
to tune 7 at a constant temperature. 

In Fig. [3] (b) and (d) , the values of Ic and Tjy result- 
ing from the fits are plotted as a function of cos(7r<I>/<I>o)i 
i.e., the normalized Ej. Ic has a linear dependence on 
cos(7r$/$o), as expected. T/v also increases monotoni- 
cally from $ = $0 to $ = 0, but the dependence of T^v 
with flux is difficult to interpret within the AH theory. 
For conventional SQUIDs, it is known that near half flux 
quanta there are multiple metastable states which may 
lead to enhanced fluctuations. However, such enhanced 
fluctuations have not been reported in previous escape 
rate experiments. I^SlEZI A somewhat similar scenario is 
the mesoscopic fluctuations in SNS junctions consider- 
ing the situation that several stationary states may exist 



at a given current 



4 38 39 



but the supercurrent fluctuation 



cannot explain the flux dependence. 

It is clear that the fitted values of T/v are larger for si 
than for s2. Since si has thinner superconducting leads, 
a larger Tjv may suggest that the thermal activation en- 
ergy barrier is related to the dimensions of the super- 
conducting leads, reminiscent of the fact that in TAPS 
theory AF is proportional to the cross-sectional area. 
Further investigation is required to verify this claim. 

In recent experimental investigations on proximity 
jimctionsjS'^ current-voltage characteristics (CVC) were 
often presented, and strong broadening was reported 
when comparing experimental data to the intrinsic 
broadening predicted by the AH theory and RSJ model. 
Here the measured CVC (measured at the same time as 
the differential resistance) are plotted in Fig. H for com- 
parison. The strong broadening can be fit well in the 
regime I < Ic hy Eq. ^ with the same Tjv in Fig. [3| At 
higher current bias the CVC collapse to a single curve, 
and cannot be fit. For sample s2, although the supercon- 
ducting electrodes are wider than that of si, the critical 
current is smaller probably due to a less transparent in- 
terface. The deviation at higher bias is less obvious for s2 



TABLE II: Fitting parameters. See Fig.[3]for fitting of dV/dl using AH theory, Fig.jljfor fitting of IVC using AH tiieory, Fig. [5] 
for fitting of IVC using LAMH theory, and Fig. p] for fitting of R{T) with three different models. 
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FIG. 4: (Color online) Current voltage characteristics at 0.4 
K for si: (a) and (b), and s2: (c) and (d). The absolute value 
of voltage is plotted as the y axis in (a) and (c), in log scale, 
and in (b) and (d), in linear scale. Solid lines are fits using 
AH theory with the same Tn in Fig. [S] 



since the spreading of Ic is smaller. We also show in Fig. 
4 the CVC on a logarithmic scale for later comparison 
with fits using the TAPS model. 

As discussed in the introduction, in terms of analy- 
sis, superconducting proximity junctions lie somewhere 
between weakly coupled conventional junctions and 
(strongly coupled) superconducting wires. Consequently, 
besides the AH theory for junctions, we also can try the 
LAMH theory for superconducting wires to fit the CVC 
of our devices. For thin superconducting wires, the aver- 
age vo ltage generated by TAPS due to thermal fluctua- 
tions idSEI 



Vl 



AMH 



— e-^^^/'^^sinh— , 
e 2kT 



(5) 



where AFq — {8y/2/3){H'^/8Tr)A^ is the energy barrier, 
ft = {L/^){AFo/kTy/^T(jl is the attempt frequency, 
tgl = 7Th/8k{Tc — T) is the Ginzburg-Landau (GL) re- 
laxation time, 5F — hi /2e is the difference in the energy 
barrier for phase slips in two directions, A is the cross- 
sectional area, and ^ is the GL coherence length. Note 
that the hyperbolic sine term in Eq. ^ has the same 
form as in the AH theory (see Eq. (4])). For our devices, 
Ai^o should be considerably lower than A_Fo for a super- 



FIG. 5: (Color online) Current voltage characteristics at 0.4 
K for si, (a), and for s2, (c), fitted by V{I) = Vo sinh(7/7o). 
In (b) and (d) the fitting parameter Vb is plotted as a function 
of fiux. 



conducting wire . Using interrelations of quantities from 
GL theory, the energy barrier can be reformulated aa^^l 



AFo = ^/&{h/2e)h, 



(6) 



which is similar to the form of Ej in Eq. (2). Thus, 
AFo/kT in the LAMH theory is comparable to 7 in the 
AH theory. 

At constant temperature, Eq. ([5]) can be simplified to 
Vlamh{I) = Vbsinh(///o), as O and AFq should have 
a much weaker dependence on the bias current / than 
the hyperbolic sine term.BH Here Vq is the flux dependent 
prefactor, and /q — ieksT/h is a constanPS at fixed 
temperature. At 0.4 K Jo is 0.0052 ^lA. In Fig. [3] CVC at 
different fields are fitted by this simplified equation (see 
Tabic [h] for fitting parameters). The fitting parameter Iq 
for si (s2) is ~ 0.07 (0.04) iiA, about 10 times larger than 
the expected value. Previously, about 40% increase of 
Iq was reported for MoGe superconducting nanowires,'^ 
but here the increase of /q is much larger. To characterize 
the increase of Iq , we can again assume an effective noise 
temperature Tjv that is about 10 times higher than the 
bath temperature T^, similar to that in fitting with AH 
theory. However, here T/v is flux independent, in contrast 
to the AH fit. 

The constant Vq can be fit with a simple exponen- 




FIG. 6; (Color online) Resistance tail at zero field for si, (a), 
and for s2, (b). The green lines are fits using Little's theory, 
the blue lines are fits using LAMH theory, and the red lines 
are fits combining Little's equation and quasiclassical theory 
for proximity junctions (RsNsiT)). The fitting parameters 
are listed in Table HIl 



tial equation y = Ae ^'^ as shown in (b) and (d) of 
Fig. [5) following Eq. ^. Using Eq. ^ and Ic{'^) = 
/c(0)]cos(7r$/$o)|j we can rewrite Vq as 



Vo{x) = Ae 



-Bx 






0), 



(7) 



where x = cos($/$o)- The factor in front of x can be 
reformulated as B = a/6/c($ = 0)/nlo. In Fig. [H] (a), if 
we define the critical current Ic as the current at which 
the measured voltage approaches 1 /zV, then at zero flux 
Ic ^ 0.56 fiA. Since Iq — 0.07 /lA for si, the expected 
B ^ 6.24, very close to the fitted value of Bfu — 6.2 
(see Table In]) . This means we do not need to introduce 
parameters other than the flux independent T^ to fit all 
CVC data. For s2 the expected B is 6.83, slightly higher 
than Bfit = 5.8, which could be due some uncertainty 
in defining the Ic values. Our analysis indicates that 
in general, it seems that the LAMH theory is better for 
fitting the CVC for proximity junctions. 



close to the Cfu extracted from the R(T), as listed in 
Table ini 

In the full LAMH theory the attempt frequency term 
is added to Eq. ^ (while the prefactor Rn is dropped), 
and the resistance below T^ derived from Eq. ([5]) in the 



zero bias limit is 



R 



Vo 



LAMH 



-AFo/kT 



- Dt-^/'^a - t)9/4exp[-c(l - tfyt], (9) 

where D is the attempt frequency constant. The total 
resistance is given by i? = {Rjj^ + Rlamh)~^ ^^ quasi- 
particles provide a parallel conduction channel to the su- 
peconducting channel. The values of c for the LAMH 
fits and and the fits to Eq. (8) are very close since R{T) 
is very sensitive to the exponential term. The attempt 
frequency constant D = {8/Tr){L/S^{0))RQy/c is formu- 
lated for a superconducting wire with length L, where 
Rq = h/Ae^ is the resistance quantum. From the fitted 
values of D and c we get i/C(0) ^ 0.01, a very small 
value as expected for proximity junctions. 

In the analysis above, the fit for R{T) is valid only in a 
narrow regime below the fitted Tc, which is due to the fact 
that A_Fo is estimated using Eq. (l6| , where for supercon- 
ducting wires near T^, Ic oc (1 — t)'^^. For SNS junctions, 
quasiclassical theory has shown that in th e long junction 
limit (A > ETh), the simplified solution ifpEI 



cRnIc = (32/(3-|-2V2))£;Th(V^T)'e-^/^^, 



(10) 



where Lt = \/hD/2nkBT is the thermal length in the 
diffusive limit, and D is the diffusion constant of the nor- 
mal metal. We can put this lc{t) dependence into AF of 
Eq. (8), which results in: 



IV. RESISTANCE TAIL 

The thermal activation theory can fit the resistance 
tail slightly below Tc for a thin superconducting wire. In 
the sim plest case, we have the Arrhenius law proposed 
by LittlePEH 



^ Rn exp[~c{l - t)^/'^ /t], 



(8) 



where t = T/Tc, c = AFf){0)/kTc is the temperature 
independent activation constant. 

As we noted earlier, the energy barrier for a prox- 
imity junction could be considerably lower than that 
for a superconducting wire, and Eq. ([6]) is used for a 
practical estimate of c from the measured Ic- Taking 
B = VGIc/ttIo ^ 6.24 obtained at 0.4 K, and the fitting 
parameter Tc = 0.78K (see Table [h] for the fitted Tc val- 
ues with Little and LAMH models, which is somewhat 
arbitrary"*^), from B = c{l — t)^^"^ jt we find c = 9.41, 



e^RN 3 + 2^/2 



Tib, 0= — , . 

^hD/2TTkBTc 



(11) 



The formula is valid to much higher temperatures as 
shown by the red lines in Fig. |6l and there is no need to 
assume an arbitrary Tc for the junction that is lower than 
Tc of Al. Close to the Tc of Al the fits deviate slightly 
from the data, probably because the long junction limit 
is no longer applicable, and also may due to the inverse 
proximity effect. With sample parameters L ^ 1 /im, and 
D ^ 100 cm^/sec, a and b are estimated to about 5 x 10** 
and 10, the same order as the fitted values. We note 
that with the other thermal activation model, similarly 
fits can be achieved by using the zero current limit of 
junction resistance in AH theory Rah = RnIo'^{i/2), 
where Xg is the modified Bessel function. However, in 
that case the fitted a values are much smaller than the 
estimate. 




FIG. 7: Schematic of the SNS SQUID with applied magnetic 
flux $ and phase difi'erence ^p between two superconducting 
reservoirs. Numbers 1-4 correspond to normal-metal wire seg- 
ments. 



V. SUMMARY 

In conclusion, we have measured nanoscale proximity 
junctions made of a mesoscopic normal-metal loop con- 
tacted by thin superconducting electrodes. The effect of 
thermal fluctuations in these devices can be characterized 
by an effective noise temperature. We find that the mea- 
surement results can be better described by the LAMH 
theory for thin superconducting wires, rather than by the 
AH theory for weak coupled Josephson junctions. With 
the LAMH theory, only the effective enhanced noise tem- 
perature (equivalent to the reduction of thermal activa- 
tion barrier) is required to fit the data at different flux 
values. We also find that quasiclassical theory can be 
combined with thermal activation theory to fit R(T) of 
proximity junction devices. These observations indicate 
that for nanoscale proximity junctions we need to con- 
sider the finite size effect of the superconducting elec- 
trodes, and a quasiclassical model more completed than 
the simple one-dimensional model presented here may be 
needed to compute the reduction of the energy barrier. 

This research was conducted with support from the 
National Science Foundation under grant No. DMR- 
1006445. 



Appendix: Simulation 

The simulations shown in Fig.l were done by solv- 
ing the lineariz ed U sadel equations in the Riccati 
parameterization.SSBSl In tlie normal-metal wire, the sim- 
plified Usadel equation reads 



(9^7 -I- 2ie7 = 0. 



(A.l) 



Here the coherent function 7 is a complex function of 
distance x and energy e, normalized by the total length L 
and Thouless energy Eth respectively. The equation can 
be readily solved with appropriate boundary conditions. 



For the SNS SQUID model shown in Fig. ^ we 
match the solutions of 4 normal-metal wires by speci- 
fying bound ary co nditions at the two nodes and at the 
SN interface j ^^ l "^^ ! e.g., at the node connecting wire 1, 2, 
and 3 we have 



7.(1) 

7.(1) 

^(1)9,7.(1) 



7/(2), 
7/(3), 
A(2)9,7;(2) 



A(3)9,7K3), (A.2) 



where A(n) is the normalized product of the cross- 
sectional area and the normal conductance of the wire, 
and the subscript ^, r indicates the left or right end of 
the wire. Similar boundary conditions are used for the 
node connecting wire 2, 3, and 4 with phase factor due 
to the flux taken into account. As mentioned in the main 
text, here we assume the gap in the superconductor is not 
affected by the normal-metal wire so that the coherent 
function takes its value in bulk superconductor. 



7o 



A 

iVA2 - £2 



(A.3) 



After solving the equations, we compute the spectral 
supercurrent 



3E = ^ 



2(797 — 7^7) 



(1+77)2 



(A.4) 



where 7 is the time reversed coherent function. Then we 
integrate the spectral supercurrent to find total super- 
current 



^s = / 3e tanh( 



2T 



)de. 



(A.5) 



Since the spectral current is a conserved quantity, this 
calculation can be done at any point of the wire. For 
wire segments 1 and 4, the total current is the super- 
current through the SNS SQUID, a dimensionless quan- 
tity normalized by EcjeB}^^^^^^ where R is the normal 
resistance of wire of length L. The critical current at 
particular $ is calculated by finding the maximum to- 
tal current while varying ip. With the leading order ap- 
proximation, the critical current shows the conventional 
|cos(7r($/<i>o))| modulation. The energy of the system as 
a function of the flux $ and phase difference ip is then 
computed according to Eq. p|. The results shown in 
Fig. 1, normalized by {h/2e^)Ec, were computed for 
following parameters: A — 30, T — 10, both normal- 
ized by Ec, length is 1/3 for all 4 segments, normalized 
by L, and A is [1,1,1,1] for the symmetric device and 
[1,1,0.75,1] for the asymmetric device, normalized by R 
and cross-section area. 
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